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Abstract 

Purpose: To investigate the correlation between the expression of Epidermal Growth Factor receptor (EGFr) and 
the reduction of the effective doubling time (7 D ) during radiotherapy treatment and also to determine the dose per 
fraction to be taken into account when the overall treatment time (OTT) is reduced in accelerated radiotherapy of 
head and neck sguamous cell carcinoma (HNSCC). 

Methods: A survey of the published papers comparing 3-years of local regional control rate (LCR) for a total of 
2162 patients treated with conventional and accelerated radiotherapy and with a pretreatment assessment of EGFr 
expression, was made. Different values of T D were obtained by a model incorporating the overall time corrected 
biologically effective dose {BED) and a 3-year clinical LCR for high and low EGFr groups of patients (H EGFr and L EGFr ), 
respectively. By obtaining the T D from the above analysis and the sub-sites' potential doubling time (T pot ) from flow 
cytometry and immunohistochemical methods, we were able to estimate the average T D for each sub-site included 
in the analysis. Moreover, the dose that would be reguired to offset the modified proliferation occurring in one day 
(Dpmiif), was estimated. 

Results: The averages of T D were 77 (27-90) 95 o /o days in L EGFr and 8.8 (7.3-1 1 .0)g 5 o /o days in H EGFr , if an onset of 
accelerated proliferation T K at day 21 was assumed. The correspondent H EGFr sub-sites' T D were 5.9 (6.6), 5.9 (6.6), 4.6 
(6.1), 14.3 (12.9) days, with respect to literature immunohistochemical (flow cytometry) data of 7" pof for Oral-Cavity, 
Oro-pharynx, Hypo-pharynx, and Larynx respectively. The D pra /,y for the H EGFr groups were 0.33 (0.29), 0.33 (0.29), 0.42 
(0.31), 0.14 (0.15) Gy/day if a = 0.3 Gy 1 and a//3 = 1 0 Gy were assumed. 

Conclusions: A higher expression of the EGFr leads to enhanced proliferation. This study allowed to quantify the 
extent of the effect which EGFr expression has in terms of reduced T D and D pro i if for each head and neck sub-site. 

Keywords: EGFr, Doubling time, Potential doubling time, Cell loss factor 



Background 

HNSCC accelerates the production of clonogenic cells 
during radiotherapy, whereby an amount of a given dose 
of radiation may be used to sterilize cells produced dur- 
ing the treatment [1]. Therefore, by maintaining the 
same total dose, a reduction of OTT results in increased 
T-site control. 
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The benefit of reduced OTT has been tested in several 
studies comparing conventional treatment with acceler- 
ated fractionation schedules. The data showed an 
improved 5-year LCR [2,3]. 

However, the response is heterogeneous with respect 
to the different expressions of EGFr in the patient popu- 
lation and also to the sub-sites, as accelerated repopula- 
tion of clonogenic tumour cells and locoregional control 
could arise. 

EGFr is overexpressed in the majority of HNSCC [4] 
and activation of the receptor leads to phosphorylation 
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of the tyrosine kinase domains on the intracellular part 
of the receptor, activating downstream cascades which 
result in altered gene activation and modulation of the 
cell products. This has been related to increased cell 
proliferation, decreased apoptotic activity, increased 
angiogenesis, increased invasive and metastatic potential, 
and hence increased resistance to anti tumour therapy. 

Furthermore, it has been shown that tumours with 
high expression of EGFr have a better LCR when treated 
with accelerated radiotherapy, while there was no benefit 
of acceleration in tumours with low EGFr. 

Consequently, high EGFr has been suggested as a 
negative prognostic factor when OTT is prolonged, and 
as a positive prognostic factor when treatment time is 
reduced [5]. 

The aim of the present study is to investigate the 
correlation between EGFr expression and the reduction 
of T D during radiotherapy treatment and also to deter- 
mine the dose per fraction to be taken into account 
when the OTT is reduced in accelerated radiotherapy of 
HNSCC. 

To achieve this goal, the data published in the litera- 
ture were reviewed and analyzed by comparing different 
3-year LCR and OTT for various dose fractionation 
schemes, also taking into account different sub-sites of 
HNSCC. 



Methods 

Literature review 

The primary end point considered for the present 
analysis was LCR, defined as the probability of avoiding 
local regional recurrence of cancer at the primary 
tumour site (T) or nodal (N) position, within 3-years 
after the end of radiotherapy. 



A survey of the published papers comparing LCR for 
patients with HNSCC treated with conventional and 
accelerated radiotherapy, respectively, and with a pretreat- 
ment assessment of EGFr expression, was made [6-11]. 

In the published papers, different criteria of EGFr 
expression assessment according to the intensity of stain- 
ing were used. EGFr expression was classified by the 
investigators, with several quantitative or semi-quantitative 
scoring systems, i.e., absent, minimal, moderate, or intense 
staining (Table 1). The main characteristics for selection 
were conventional and accelerated fractionations, different 
OTT, assessment of EGFr expression and LCR, as listed 
in Table 1. 

Only those studies which reported a median follow-up 
of at least 3-years were included in the analysis. Table 2 
lists the main clinical characteristics of the patients, 
namely age, sex, primary site, T stage and N stage. Fur- 
ther clinical information are in the reviewed papers. 

Radiobiological analysis 

The tumour effects were evaluated by the overall time 
corrected BED as in eq. (1) 



BED = nd 1 



_d_ 



ln 2 ( r -^) 



(1) 



where n is the number of fractions of size d in Gy, a and 
/? are the linear quadratic coefficients of dose, T is the 
overall time, is the onset time for accelerated prolif- 
eration and T D the effective doubling time. The first 
term in eq. (1) (the dosimetric component, see 
Appendix A), is affected by differences in EGFr expres- 
sion because of modification to a and 8 parameters that 
describe the intrinsic and repair radiosensitivity of 
tumour types, respectively. We add the subscripts H and 



Table 1 Treatment characteristics of the selected group 


Author 


Fractionation 


OTT(days) 


SI H % 


SI L % 


LCR H 


LCR L 


HR 


S 


Er ksen AO a [6] 


33x2Gy 


66 


SI > 50 


SI < 50 


0.15 


0.44 


Y 


N 


Er ksen AO b [6] 


33x2Gy 


45 


SI > 50 


SI -c 50 


0.64 


0.55 


Y 


N 


Er ksen AO c [6] 


33x2Gy 


38 


SI > 50 


SI -c 50 


0.77 


0.57 


Y 


N 


Er ksen RO a [7] 


33x2Gy 


45 


SI > 50 


SI < 50 


0.57 


0.63 


Y 


N 


Er ksen RO b [7] 


33x2Gy 


38 


SI > 50 


SI < 50 


0.70 


0.62 


Y 


N 


Bentzen JCO a [8] 


33x2Gy 


45 


SIS 40 


SI < 40 


0.30 


0.45 


N 


N 


Bentzen JCO b [8] 


36x1.5 Gy 


12 


SI >40 


SI < 40 


0.54 


0.49 


N 


N 


Suwinski IJROBP a [9] 


35x1.8 Gy 


47 


SI > 33 


SI < 33 


0.33 


0.70 


N 


Y 


Suwinski IJROBP b [9] 


35x1.8 Gy 


35 


SI > 33 


SI < 33 


0.58 


0.73 


N 


Y 


Smid IJROBP a [10] 


25x2Gy + 5x2.5 Gy 


~ 46 


A/m 


M/l 


0.69 


0.65 


N 


Y 


Smid IJROBP b [10] 


25x2Gy + 5x2.5 Gy 


~ 34 


A/m 


M/l 


0.91 


0.68 


N 


Y 


Chung IJROBP a [11] 


35x2Gy 


4/ 


SI > 80 


SI < 80 


0.36 


0.61 


N 


N 


Chung IJROBP b [11] 


30x1.8 Gy+ 12x1.5 Gy 


38 


SI > 80 


SI < 80 


0.54 


0.68 


N 


N 



Abbreviations: S/ = Staining Intensity cut-point; A/m = Absent/minimal; M/l = Moderate/Intense; HR = Hypoxic Radiosensitizer; S=Surgery. 
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Table 2 Clinical characteristics of the selected groups 



Author 




E1 


E2 


BE 


su 


SM 


CH 


number of patients 


209 


803 


304 


148 


165 


533 


Sex % 


M 


74.2 


81.0 


72.0 


90.0 


66.0 


79.2 




F 


25.8 


19.0 


28.0 


10.0 


34.0 


20.8 


Primary site % 


Oral Cavity 


/ 


12.0 


13.0 


50.6 


100.0 


10.3 




Oropharynx 


/ 


52.5 


28.0 




/ 


60.4 




Hypopharynx 


/ 




12.0 


/ 


/ 


13.7 




Nasopharynx 


/ 


/ 


3.0 


/ 


/ 


/ 




Larynx 


100.0 


35.5 


44.0 


49.4 


/ 


15.6 


T stage % 


T1 


6.0 


67.0 


3.0 


26.1 


2.0 


5.8 




T2 


37.0 




42.0 




25.0 


27.6 




T3 


35.0 


33.0 


33.0 


73.9 


28.0 


37.1 




T4 


11 n 
zz.u 




ii n 
zz.u 




jj.U 


1Q 1 




Tx 


/ 


/ 


/ 


/ 


/ 


0.2 


N stage % 


NO 


65.0 


65.0 


63.0 


33.8 


26.0 


21.8 




N1 


35.0 


35.0 


18.0 


66.2 


16.0 


18.6 




N2a 






15.0 




1.0 


9.6 




N2b 










39.0 


17.8 




N2c 










17.0 


19.3 




N3 






4.0 




1.0 


12.9 



Abbreviations: Ei = Eriksen-AO, E2 = Eriksen RO, BE = Bentzen JCO, SU = Suwinski 
IJROBP, S/M = Smid IJROBP, CH = Chung IJROBP. 



L to indicate high or low EGFr expression respectively 
(BED H (d) or BED L (d)). The second term (the temporal 
component, see Appendix A) is affected by differences 
in EGFr expression due to the presence of the a param- 
eter [pc H or a L ) and T D (T DH or T DL ). Superscripts S and 
F are specified to distinguish between conventional 
(S = Slow) and accelerated (F=Fast) fractionations, 
respectively. 

From BED we have the standard model of tumour 
control probability (TCP) using the linear-quadratic 
model incorporating the Poisson's low [12], 

TCP= exp(-N-S) 

= exp[-7V ■ exp(-a ■ BED)] (2) 

where N=p-V (/3 = cell density and V= volume) repre- 
sents the initial number of potential proliferating cells in 
the tumour. Therefore, the cell survival probability being 
S = exp(-a-BED), the TCP represents the probability of 
avoiding local recurrence [13] at total dose D = n-d 
whereby we write TCP = LCR. 

Moreoever, in order to analyze the effects of EGFr ex- 
pression due to the change in the OTT, the papers 
chosen in the survey had the same dose per fraction and 
total dose but a different OTT. 



Thus, by taking the natural logarithms of eq. (2) writ- 
ten for fast and slow fractionations, dividing the result- 
ant equations and by taking the natural logarithm again, 
we get 




for low EGFr expression group (see Appendix A). 

This expedient allows to eliminates the dependence of 
findings from the choice of dose fractionation and from 
the estimated values of a and B. The equations (3.a) and 
(3.b) are also independent notwithstanding the assump- 
tion about number of cells N. The uncertainties arising 
from these assumptions strongly influences the results 
of the other models that depend on such parameters. 
Therefore, this is the main advantage of equations (3.a) 
and (3.b). 

In each of these equations appears only one unknown 
(the effective doubling time) for which, being in a linear 
form, they are suitable for an easy comparison between 
LCR due to different EGFr expression groups with dif- 
ferent OTT. This assessment was done by evaluating the 
differences of angular coefficients (\r\2/T DH vs \r\2/T DL ) 
from the correspondent regression lines obtained by 
LCR available in literature (Figure 1). For those papers, 
where in addition to differences of OTT there are also 
differences in terms of dose fractionation, the correction 
as described in Appendix B was done. 

Furthermore, dividing equations (3. a) and (3.b), we 
also obtained the ratio of the actual doubling times be- 
tween the H EGFr and L EGFr groups that allows a direct 
analysis of the EGFr effects (Figure 2) as follows 

Clinical analysis 

The actual doubling times obtained from the above ana- 
lysis, represent a weighted average of the doubling times 
from different sub-sites as oral cavity (18.2% of patients), 
oro-pharynx (30.3% of patients), hypo-pharynx (14.8% of 
patients), and larynx (36.4% of patients). 

These sub-sites contribute differently to the average 
T D because they have different T pot . However, the T D 
for each sub-sites can be estimated if T pot and the cell 
loss factor (tp) are known as described by Steel [14]. In 
particular T pot can be measured by a single biopsy 
with flow cytometry as well as immunohistochemistry 
techniques. 
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Figure 1 Linear regressions for L EGFr group (diamonds) and H EGFr group (circles) with y = /nf/nLCR s /lnLCR F j and x = (T s -T f ) (days). The 

higher angular coefficient of H EGFr group line demonstrate a reduction of effective doubling time with respect to L EGFr group (proportionality 
with the inverse of the effective doubling time). Sub-figures refer to a different onset of accelerated repopulation (a) 1^ = 0, (b) 7^= 14, (c) 7^=21 
(d) T k = 28 days. Error bars represent the 95% confidence intervals. 



Therefore, the average cell loss factor was estimated 
using pretreatment data about T pot available in literature 
[15,16], then the actual doubling time for each i sub-site 
{'T D ) was obtained (see Appendix C). 

Moreover, from l T D we also estimated the dose (in 
fractions of size d) that would be required to offset the 
effect of proliferation occurring in one day [17] by the 
follows equation 



•D 



prolif 



ln2 

T T D (a + ]ki) 



(5) 



Statistical analysis 

In all the original studies of the survey the primary end- 
point was LCR, 3 or 5-years after completion of radio- 
therapy, although only the 3-year LCR were extrapolated 
in order to compare the homogeneous parameters. LCR 
were assessed by the Kaplan-Meier method with a log 
rank test (statistical significance: p < 0.05, two-sided). 
The LCR 95% confidence intervals are obtained by 
Greenwood' formula [18]. Comparison between regres- 
sion lines was done by Fisher's exact test. 



Results 

Table 1 and 2 describe the main clinical characteristics 
and treatment parameters of the selected groups in the 
survey. 

Linear regression lines from equations (3.a) and (3.b), for 
H EGFr and L EGFr groups, are shown in Figure 1 with respect 
to different choices of the onset for accelerated repopula- 
tion (2a). The significant distinction of the angular coeffi- 
cients for different groups (p-values < 0.02) correspond to 
an average T D of 77 days (27-90) 95% for L EGFr and to an 
average of 8.8 days (7.3-11.0) 95% for H EGFn if an onset of 
accelerated proliferation T K at day 21 was assumed. 

In Figure 2 the significant H EGFr T D reduction with re- 
spect to L EGFr T D for each head and neck sub-site, are 
shown by varying T K . 

In Figure 3 the averages of D pro uf are shown based on 
the flow cytometry and immunohistochemical methods 
to estimate the sub-sites 'T pot . The maximum value of 
Dproitfi up to about 0.5 Gy/day, is obtained corresponding 
to an onset of accelerated repopulation that starts from 
the fourth week {T K at about 28 th day). Sensitivity ana- 
lysis is shown with respect to different values of a with 
a/fi = 10 Gy. 

The weighted average potential doubling times < l Tp 0 t> 
of 5.2 days and 3.4 days [15,16] were obtained 
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corresponding to averages for cell loss factors as < 
0(fcm)>= 0.41 (0.29-0.52) 95 o /o and<0 (H i)>= 0.61 (0.53- 
0.69) 95 o /o with respect to the flow cytometry and immu- 
nohistochemistry, respectively. 

Table 3 reports numerical results for each sub-site in 
H EGFr group with 'T D and 'D proU f calculated for different 
values of a (a//? = 10 Gy). It may be noted that the l T D 
for each i sub-site is almost twice of l T pot obtained by 
flow cytometry and more than double of l T pot obtained 
by immunohistochemistry. This means that a pre- 
treatment assessment of D pro uf by flow cytometry or 
immunohistochemistry may significantly overestimate 
the dose required to offset the accelerated proliferation 
occurring in one day. 

In Figure 4, the histogram of the ratio between T DL 
and T DH (eq. 4) shows an average reduction of about 7 
times in average (6.6-8.3) 95% for the H EGFr group with 
respect to the L EGFr . This ratio could have significant 
implications on the clinical management of these patient 
groups. In fact, while the H EGFr group would benefit from 
an increase of the dose/fraction (Hypo-fractionation) and 
the consequent reduction of OTT to compensate for 
the increase in the proliferation rate - corresponding to 
a reduced T D -, the L EGFr group does not require a 
reduction of OTT for which it would be more indicated a 
reduction of the dose/fraction (Hyper- fractionation) which 
would result in a reduced toxicity for all the organs at risk. 



Discussion 

In the recent years there has been a great interest to find 
factors that predict tumours suitable for accelerated 
radiotherapy and considerable interest has been given to 
cell kinetic parameters such as the T pot . Since regener- 
ation and tumour cell proliferation are mechanisms at 
the cellular level, particular attention has been focused 
on identifying the specific cellular characteristics, such 
as variations in EGFr expression. The latter is an import- 
ant mediator of cell growth and its over-expression has 
been associated with tumour progression and poor sur- 
vival in many solid cancers. Several studies have demon- 
strate the potential of EGFr as a predictive and 
prognostic marker in radiotherapy for HNSCC [19]. 

In the present study a direct demonstration of the link 
between EGFr status and the time factor in fractionated 
radiotherapy, has been made. All the clinical studies sur- 
veyed, from the available literature, had a random allocation 
for "reduced" or "conventional" OTT and demonstrated an 
increase in LCR when the OTT was reduced. 

Unfortunately, OTT reduction yields clinical benefits 
in terms of LCR but could worsen the radiation-induced 
acute side effects which need to be carefully evaluated 
using appropriate radiobiological models [20]. 

Moreover, some studies also demonstrated that 
tumours with high EGFr respond better to the reduction 
of the OTT compared to low EGFr tumours [6,7,9]. The 



Table 3 Numerical results of sub-site 



Sub-site Number Doubling Flow Cytometric (FCM) Immunohistochemical (Hi) 





of 
patients 


time 


days 


Dp (0.25) 


Dp (0.3) 


Dp (0.35) 


days 


Dp (0.25) 


Dp (0.3) 


Dp (0.35) 


Oral Cavity 


394 (18.2%) 


T D 


6.6 (5.5-8.2) 


0.35 (0.28-0.42) 


0.29 (0.23-0.35) 


0.25 (0.20-0.30) 


5.9 (4.9-7.4) 


0.39 (0.31-0.47) 


0.33 (0.26-0.39) 


0.28 (0.22-0.34) 






Tpot 


3.9 


0.59 


0.49 


0.42 


2.3 


1.00 


0.84 


0.72 


Oropharynx 


655 (30.3%) 


T D 


6.6 (5.5-8.2) 


0.35 (0.28-0.42) 


0.29 (0.23-0.35) 


0.25 (0.20-0.30) 


5.9 (4.9-7.4) 


0.39 (0.31-0.47) 


0.33 (0.26-0.39) 


0.28 (0.22-0.34) 






Tpot 


3.9 


0.59 


0.49 


0.42 


2.3 


1.00 


0.84 


0.72 


Hypopharynx 


320 (14.8%) 


T D 


6.1 (5.1-7.6) 


0.38 (0.30-0.45) 


0.31 (0.25-0.38) 


0.27 (0.22-0.32) 


4.6 (3.9-5.8) 


0.50 (0.40-0.59) 


0.42 (0.33-0.49) 


0.36 (0.28-0.42) 






Tpot 


3.6 


0.64 


0.53 


0.46 


1.8 


1.28 


1.07 


0.92 


Larynx 


784 (36.3%) 


Td 


12.9 (10.7-16.0) 


0.18 (0.14-0.22) 


0.15 (0.12-0.18) 


0.13 (0.10-0.15) 


14.3 (12.4-18.1) 


0.16 (0.13-0.19) 


0.14 (0.11-0.16) 


0.12 (0.9-0.13) 






Tpot 


7.6 


0.30 


0.25 


0.46 


5.6 


0.41 


0.34 


0.29 



T D (days) and D protj f (Gy/day) for high EGFr expression group (<'T pot >=5.2 days by flow cytometry; <T pot > =3.4 days by immunohistochemistry; T K = 2'\ days). The 'T D for each / sub-site (upper rows) is almost twice of 
'T pot by flow cytometry and more than double of 'T pot by immunohistochemistry from literature (bottom rows) [1 5]. Consequently, a preclinical estimation of 0 O{O y lf by flow cytometry (left columns, bottom rows) or 
immunohistochemistry (right columns, bottom rows), results almost double and more than double if compared to the estimation of D praW obtained by mathematical model (left and right columns, upper rows), 
respectively. The 95 % confidence intervals are shown within the brackets. 

Abbreviations: 'D pt0 2 s ; 0.3,- 0.35) = Opra/tf calculated with a = 0.25, a = 0.3, a = 0.35 Gy 1 respectively, a/fi = 1 0 Gy, d = 2 Gy. 
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Figure 3 Curves of D prom obtained by varying T k , for H EGFr (continuous) and L EGFr (dashed) groups based on the flow cytometry (left 
column) and immunohistochemistry (right column) methods to estimate the sub-sites 'T pot . Sensitivity analysis is shown with respect to 
different values of a with a/j8= 10 Gy: a = 0.25 Gy" 1 (upper curves); a = 0.3 Gy" 1 (mean curves); a = 0.35 Gy" 1 (bottom curves). Differences are shown 
for Oral-Cavity (a,b), Oropharynx (c,d), Hypopharynx (e,f) and Larynx (g,h). 
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CH 



Figure 4 Histogram of ratios between T D values in the l_EGFr and HEGFr groups. The error bars represent the 95% confidence intervals while 
the overlay trend line represent the average of values. Abbreviations: E1 =Eriksen from AO, E2 = Eriksen RO, BE = Bentzen JCO, SU = Suwinski 
IJROBP, SM = Smid IJROBP, CH = Chung IJROBP. 



response was heterogeneous if referring to the sub-sites 
included in the analysis. 

Therefore, our intent was to evaluate the extent 
of accelerated proliferation due to an EGFr over- 
expression, in terms of reduced actual doubling time as 
well as required dose to offset the effect of proliferation 
occurring in one day. 

To obtain these results no assumption was made with 
the exception of the validity of a linear quadratic and 
TCP model. Therefore, the fact that the EGFr expression 
changes the radiosensitivity and the proliferation rate of 
the cells, has to be necessarily included in these models 
as a variation of the parameters {a, /S and T D ) describing 
them. 

Although in all the studies selected for the survey, the 
accelerated repopulation of tumour cells during radiother- 
apy was suggested as an important cause of treatment 
failure, the main difficulty in our analysis was that head 
and neck cancer represents a heterogeneous group of 
cancers and the benefit is not act equal for the different 
tumours. We also attempted to estimate these differences. 

We are aware that an important drawback in the ana- 
lysis is to be found the differences among treatment mo- 
dalities. Some studies have included radiotherapy alone, 
others postoperative radiotherapy, others the use of 
radio-sensitizing hypoxic drugs. We therefore stress the 
versatility and enormous potential of the method we 
propose. 

Indeed, the relationship between elements representa- 
tive of the radiation effects calculated only on groups of 
patients who undergo the same treatment, is based 



accordingly. In other words, although radiotherapy alone 
is profoundly different from postoperative radiotherapy 
or from radiotherapy combined with radio-sensitizing 
drugs, the relationship of the effects calculated within 
the same type of therapy, nullify these contributions, 
allowing to obtain only those due to the different expres- 
sion of EGFr. 

The validity of this statement is confirmed by the very 
low dispersion of data around the linear regression lines 
obtained from them. Our results clearly demonstrate 
proportionality between differences in treatment dur- 
ation and correspondent ratios of LCR (the latter in 
logarithmic form). This strong linearity allowed us to 
quantify the reduction of actual doubling time of the 
H EGFr group with respect to the L EGFr group. 

Unfortunately, in the papers, different definitions for 
the level of cut-points of EGFr expression were used. In 
some studies a cut-point of 50% was chosen as being ob- 
jective and reproducible, others fixed a cut-point of 33%, 
40%, 80%, etc. However, it is obvious that no dichotom- 
ous division between high and low EGFr tumour expres- 
sion exist and a continuous variable must apply. In 
addition, because samples for the various studies were 
collected from various pathology departments and stain- 
ing intensity can be dependent on tissue fixation [21], 
the evaluation of staining intensity was not entirely 
homogeneous. This was certainly the greatest source of 
approximation in the quantitative results obtained. 

Despite these limitations, our results indicate a clear 
reduction of effective doubling time T D in H EGFr with re- 
spect to the L EGFr groups. This reduction did not so 
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excessive as the necessary to reach the minimum value 
which is represented by T pot (that is the limit where the 
cell loss fraction is reduced to zero and the proliferation 
is fastest). 

We consider this result very important, especially be- 
cause an accurate estimate of T D allows to obtain the 
equivalent dose for accelerated repopulation that is es- 
sential to making rational adjustments to the overall 
dose when the overall time is increased. This has be- 
come more than just an academic question in the area 
of IMRT when, instead of using shrinking field techni- 
ques, radiation oncologists commonly use a differential 
dose per fraction to deliver graded doses in the same 
overall treatment time. 

Our results consisted of D pro uf systematically lower 
than those accepted in the literature that are often 
obtained through an evaluation of-the potential doubling 
time, which is a characteristic of each proliferative cell, 
and not through the effective doubling time, which is-a 
characteristic of a group of cells [22] . 

In the case of oropharynx, for example, we obtained 
values up to 0.39 Gy/day (0.31-0.47) 95% while in the lit- 
erature we have values between 0.48-0.68 Gy/day 
[23,24]. For the larynx we obtained values up to 
0.16 Gy/day (0.13-0.19)95%, while other estimates for this 
tumour are between 0.3 and 0.5 Gy/day [17]. Only in the 
case of hypopharynx we had values greater than 0.5 Gy/ 
day (0.40-0.59) 95 %- 

The difference can partly be explained by the hetero- 
geneous behavior of the different sub-sites involved in 
our analysis (a specific sub-site clinical study could dis- 
criminate more finely between different contributions). 
However, our opinion is mainly based on the interpret- 
ation of the correlation between T D and T pot . 

As a first hypothesis, the reduction of T D can be easily 
explained with a correspondent reduction of T pot , but 
clinical data has shown that for patients with short T pot 
(fast tumours) there was no statistically significant trend 
to do worse [25]. Moreover, we found a reduction of T D 
with an average factor of about 7 in the H EGFr with re- 
spect to the L EGFr group, and the same extent was never 
found by measures that assess T pot from biopsy among 
patients. 

Consequently, given that a shorter T D may also result 
from a reduced cp after the beginning of treatment (see 
Appendix C), our results suggest that the latter possibil- 
ity is favored. In this case, the tendency of cp toward 
zero, indicates a reduction of the clonogen doubling 
time T D until it equals the pretreatment T pot . Hence, our 
results for T D can be easily explained from an incom- 
plete reduction of <p toward zero. 

Furthermore, a <p reduction being associated with a 
low differentiation, would correspond to an increase in a 
non-differentiated component. 



Thus, the question arises about how two different 
results may be reconciled. 

On one hand, the simultaneous expression of a differ- 
entiated pattern and high levels of EGFr display a higher 
degree of accelerated repopulation compared to carcin- 
omas with low levels of EGFr or poor differentiation [5]. 
On the other hand, as is clear from the Steel's formula, 
the reduction of T D is due to a reduction of the differen- 
tiation levels. 

A possible explanation could be that two different 
levels of differentiation may coexist locally. 

This hypothesis is based on the clinical observation 
that high levels of EGFr expression were found to be 
more pronounced at the tumour borders compared to 
the central parts of the tumour tissue (p < 0.0001) [6]. 
Therefore, on the border of the tumour, the EGFr over- 
expression would be compatible with a low level of dif- 
ferentiation and rapid tumour growth (as from Steel's 
formula). In more central tumour areas, the low EGFr 
expression may be compatible with a high level of differ- 
entiation and reduced tumour cell proliferation. This 
spatial non-uniformity, suggests that the precise location 
of biopsy sampling and a subsequent classification of 
tumours (high or low EGFr and level of differentiation) 
are crucial. A such hypothesis, of course, requires fur- 
ther investigation in clinical studies. 

Conclusion 

Increased expression of the EGFr can lead to enhanced 
proliferation which can be countervailed by reducing the 
time available for tumour cell proliferation, thereby re- 
ducing the overall treatment time. In this case, the im- 
pact of high EGFr expression changes from being a 
negative to a positive prognostic value in terms of local 
control rate. 

In this study we introduced a model that allows to 
quantify the influence of EGFr expression in terms of 
reduced doubling time during the treatment and also 
the dose per fraction to be taken into account when 
the overall treatment time is reduced in accelerated 
radiotherapy. Furthermore, using this model, we can 
also estimate the parameters inherent in different 
sub-sites which may identify the optimal dose frac- 
tionation regime more likely to benefit these sub-sets 
of patients. 

Appendix A 

To simplify the radiobiological analysis, eq. (1) can be 
rewritten by considering the BED as the difference be- 
tween a dosimetric and a temporal component: 



BED = BED(d) - BED(T) 
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The BED(d) for high EGFr expression group, for in- 
stance, is [26]. 



BED H (d) =nd[\ 



while the BED(T) for the same group is 



BED H (T) 



ln2 



a, i 



1 d 



Yh_ 
cch 



(T-T K 



Thus, to take into account the differences of radiosen- 
sitivity as well as OTT, two and four possible expres- 
sions of BED(d) and BED(T) are considered, respectively: 

BED H (d),BED L (d),BEDH(T s ) ! BED H (T F ), 
BED L (T S ),BED L (T F ) 

For different EGFr expressions and OTT we have also 

LCR F H , LCR F H ,LCRl,LCRl 

where, for example 

exp[- P Vexp(a H (BED H (d) - BED„(T S )))\ 



LCR S H 



Therefore, by taking the natural logarithms of this ex- 
pression and dividing it for the same expression with a 
different OTT, we can nullify the contributions of pV 
and BED(d) - because of the same fractionation - 
obtaining only those due to the different expression of 
EGFr with respect to the OTT. 

For the high EGFr expression group, for instance, we 
have 

InLCRfj I \nLCR F H = exp | a H (BED H ( T s ) 
-BED H (T F ))\ 

from which, by taking the natural logarithms again, we 
have eq. (3. a). The same procedure leads to eq. (3.b). 

Appendix B 

The equations (3. a) and (3.b) are valid if the hypothesis 
of equal dosimetric component of BED in conventional 
and accelerated fractionation, is valid. However, this is 
not true for all the papers in the literature surveyed. The 
differences in terms of dose fractionation were corrected 
using the follows 

\nLCR c = e -A*B*W-*>&W) \nLCR„ c 

where the indexes nc and c stand for "non-corrected" 
and "corrected", respectively. BED s (d) and BED F (d) refer 
to different BED dosimetric components in conventional 
and accelerated fractionation, respectively. The exponen- 
tial factor incorporates the difference of BED only due 



to the dosimetric BED component, and therefore enables 
the contribution to be corrected, thanks to this 
component. 

Appendix C 

In order to estimate the actual doubling time for each 
sub-site in the analysis, different potential doubling 
times T pot were considered from literature [15]. 

The latter has been introduced by Steel as the clono- 
gen doubling time that would be measured if cell loss 
was ignored, i.e. if both daughter cells remained clono- 
genic after mitosis [14]. 

In practice clonogens are lost through many possible 
mechanisms, including differentiation, death, and metas- 
tasis, and the net result is that T D will be longer than 
T 

± pot- 

Steel's formula can be written as follows: 



1 pot 



wherein 0 is the cell loss factor. This equation shows 
that T D can be calculated if T pot and (p are known. 

In particular T pot can be measured by a single biopsy 
using flow cytometry or immunohistochemistry techni- 
ques, while the average cell loss factor < cp > was obtained 
in our analysis taking the average potential doubling 
time weighted on percentages for any sub-site (<'T pot > = 
Yli Pr'T pot with i = l,..,4 and pj = 21% for oral cavity, 
P2 = 20% for oro-pharynx, P3 = 17% for hypo-pharynx 
and p 4 = 42% for larynx), by the follows 



< 0 >= 1 



<'T pot > 
<>T n > 



The actual doubling time for each i sub-site, was then 
obtained as follows 



± pot 
< 0 > 



Results were reported in Table 3. 
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